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Abstract 

The Roy equation in the single channel case is a nonlinear, singular inte- 
gral equation for the phase shift in the low-energy region. We first investi- 
gate the infinitesimal neighborhood of a given solution, and then present 
explicit expressions for amplitudes that satisfy the nonlinear equation 
exactly. These amplitudes contain free parameters that render the non- 
uniqueness of the solution manifest. They display, however, an unphysical 
singularity at the upper end of the interval considered. This singularity 
disappears and uniqueness is achieved if one uses analyticity properties 
of the amplitudes that are not encoded in the Roy equation. 
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1 Introduction 



The elastic itit amplitude has recently been evaluated in the framework of chiral 
perturbation theory Jl| to two loops 0. The representation is valid in the low- 
energy region, where the centre-of-mass energy of the pions is less than about 
400 MeV. On the other hand, precise experimental data are presently only 
available above ~ 600 MeV. In order to connect the two regions, one may rely 
on a set of dispersion relations for the partial wave amplitudes due to Roy 0. 
These allow one to extrapolate the data down to threshold and to merge with 
the chiral expansion [[| . In the present article, we investigate this extrapolation 
procedure in the one-channel case. 

Roy's representation || for the partial wave amplitudes t\ of ttit scattering 
reads as follows, 

2 oo />oo ; ; 

Retj(s) = fcf(s) + £ E f dxKtf (s,x) Imtf, (x) ; 4 < s < 60 , (1.1) 

where / and I denote isospin and angular momentum, respectively^]. The linear 
subtraction polynomials kj are expressed in terms of the two S-wave scatter- 
ing lengths. The kernels K\y contain a diagonal, singular Cauchy kernel that 
generates the right-hand cut in the partial wave amplitudes, as well as a loga- 
rithmically singular piece that generates the left-hand cut. 

The relations ( |Llp are consequences of the exact analyticity properties of 
the 7T7T scattering amplitude, of the Froissart bound and of crossing symmetry. 
They demonstrate that the scattering amplitudes are fully determined by the 
imaginary parts of the partial waves, except for the two scattering lengths that 
play the role of subtraction constants. Combined with unitarity, ( |1.1|) amounts 



to an infinite system of coupled, singular integral equations - the Roy equations 
- for the phase shifts in a low-energy interval from threshold s = 4 to a matching 
point So < 60. In this framework, the phase shifts above the matching point Sq, 
the absorption parameters and the two S-wave scattering lengths are assumed 
to be externally assigned. The mathematical problem consists in solving the 
Roy equations with this input. 

Soon after the original article of Roy M had appeared, extensive phenomeno- 
logical applications were performed ||, resulting in a detailed analysis and ex- 
ploitation of the then available experimental data on tttc scattering. For a recent 
review of those results, we refer the reader to the article by Morgan and Pen- 
nington || . Parallel to these phenomenological applications, the very structure 
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of the Roy equations was investigated. In || , extensions of the equations ( [Lip 
were presented, valid in the larger range —28 < s < 125.31. Further, the mani- 
fold of solutions of Roy's equations was investigated as well, both in the single 
channel |5], [ll], |11[ as well as in the coupled channel case [ I2|] . In the late seven- 



ties, Pool [13|] provided a proof that the original, infinite set of integral equations 



does have at least one solution for s < 23.31, provided that the driving terms 
are not too large, see also |14[]. Heemskerk and Pool also examined numerically 
the solutions of the Roy equations, both by solving the N equation and by 



using an iterative method [15 |. 



It emerged from these investigations that - for a given input of absorptive 
parts, absorption parameters and S-wave scattering lengths - there are in gen- 
eral many possible solutions to the Roy equations. This non-uniqueness is due 
to the singular Cauchy kernel in the right-hand side of ( |1.1|) . In order to inves- 
tigate the uniqueness properties of the Roy system, one may - in a first step - 
keep only this part of the kernels, as a result of which the integral equations de- 
couple: one is left with a single channel problem, i.e. a single partial wave, that 
has, moreover, no left-hand cut. This mathematical problem was examined by 
Pomponiu and Wanders ||. Investigating the infinitesimal neighborhood of a 
given solution, they found that the multiplicity of the solution increases by one 
whenever the value of the phase shift at the matching point goes through a mul- 
tiple of 7r/2. By contrast, the number of parameters in the usual partial wave 
equation increases in general by two whenever the phase shift at infinity passes 
through a positive integer of tt, see e.g. [|H| 17]] and references cited therein. 



Atkinson and Warnock [|T0[ investigated a one-channel problem by using 
N/D methods, that do not require linearization. Their work may be sum- 
marized as follows. First, they find that solutions of the Roy equations with 
${so) > — f are members of solution manifolds depending on m parameters, 
where m denotes the integer part of 2<5(s )/vr. Second, these solutions may 
be computed through an integral equation in which the m parameters appear 
explicitly. An exception to these statements - whose proof is rather complex 
- could occur if a certain Fredholm operator had unit eigenvalue. Although 
the relevant integral equations in |Tl| exhibit the m parameters, it had not 
been shown that these are really effective: an arbitrarily chosen set may lead 
to ghosts. 

After 1980, interest in the Roy equations waned. In recent years, it has 
however become clear that there are good reasons to revive these techniques. 
First, new Ku experiments are planned [0] or are already under way |Tj| |19| . 
These will provide new information on the low-energy tttt scattering amplitude. 
A reliable analysis of available low-energy data and high-energy phase shifts 
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should then again be based on Roy equations. Second, one can establish at 



this point contact with CHPT, that will allow one [p], |21| to gain insight into 
the structure of the underlying theory, the Standard Model. First steps in 
this program have already been performed p2, p3"| . Of course, this enterprise 



requires that one understands the structure of the Roy equations in all details. 
It is the aim of this and a following article [21] to fill existing gaps in this respect 



and to provide insights into the problem from an actual point of view, see also 
Ref. 



We concentrate in this work on the uniqueness properties of the solutions 
in the one-channel case, and on their singularity structure. We start by an- 
alyzing the Roy equation with the linearization method proposed in || and 
show how the emerging integral equation can be solved by transforming it into 
a homogeneous Hilbert problem. This method is very efficient in determining 
the multiplicity of solutions. It does, however, assume the existence of a solu- 
tion. Furthermore, it is not easy to show that each solution of the linearized 
equation approximates a solution of the nonlinear one. For these reasons, we 
also investigate the Roy equation for a special class of inputs, that allow us to 
find explicit solutions. These solutions do contain parameters that exhibit their 
non-uniqueness. We then investigate the role of the unphysical singularity that 
shows up in the solutions. This singularity manifests itself as a cusp in the 
real and imaginary parts of the amplitude at the matching point. By making 
use of analyticity properties of the amplitude that are not explicit in the Roy 
equation, we find that there is a unique solution, devoid of cusps. 

Our article is organized as follows. In section 2, we formulate the math- 
ematical problem we are concerned with here. In section 3, we analyze the 
infinitesimal neighborhood of a solution. The construction of a second exact 
solution - once a first solution is known - is reduced to a linear problem in 
section 4. This allows us to give in section 5 a class of explicit exact solutions. 
In section 6 we show how uniqueness is obtained by use of analyticity proper- 
ties of the amplitudes that are not encoded in the Roy equation. A summary 
and concluding remarks are given in section 7. Details on precise mathematical 
formulations are relegated to appendix A, where we also construct the general 
solution of the linearized problem of section 2. Appendix B provides the con- 
nection with the N/D method and appendix C contains the proof of the 
uniqueness property for an analytic input. 

In order to help orient the reader, we note that the key statements in this 
article are summarized in propositions 1, 2, 3 and 4 in sections 3, 4 and 6. 
Propositions 1 and 3 were established long ago in Refs. and , respectively, 



whereas propositions 2 and 4 are - to the best of our knowledge - new results. 
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2 The one— channel Roy equation 



In order to study the non-uniqueness properties of the Roy equation, we keep 
only the diagonal, singular Cauchy kernel in (|1.1| ). The partial wave relations 
then decouple, and the left-hand cut in the amplitudes disappears. We therefore 
first explore the set of complex amplitudes / : [4, oo) — > C with the following 
properties: 

i) In an interval [4, sq] containing the threshold s = 4 and a matching point sq, 
the real part is given by a dispersion relation 

,lr°°dxlm f(x) 
Re / S ) = a + (a-4 -i — . 2.1 

7TJ4 X — 4 X — S 

ii) The imaginary part Im/ is a given input function A above so, 

Im/(s) = A(s), s>s . (2.2) 

iii) Elastic unitarity holds below sq, 

lmf(s) = a(s)\f(s)\ 2 , se[4,s ] ; a(s) = [1 - 4/s] 1 / 2 . (2.3) 

We relegate a precise formulation of the regularity properties of / to appendix 
A and simply note that, as a minimal requirement, the imaginary part Im/ 
must be continuous in [4, s ], in particular, 

lim Im f(s) = A(s ) . (2.4) 

s/s 



Equations ( (2.1| )-( fO|) constitute the mathematical problem that we discuss in 
the following: determine the amplitudes / which verify these equations for given 
scattering length a and given absorptive part A. This one-channel problem 
allows detailed analytical and numerical calculations and provides useful insight 
into the solutions of the coupled system considered e.g. in Refs. [|, |^, |T2|| . 

Elastic unitarity allows one to reduce the problem to the determination of 
a single real function on the real interval [4, s ] , because / is parametrized by 
its real phase shift 5 below s : 

/( s ) = -Le w W sin 6(s), se[A,s }. (2.5) 
a[s) 

We choose the normalization of 5 such that it vanishes at threshold, 5(4) = 0. 
The boundary condition ( |2.4j ) becomes 

sm 2 5(s ) = a(s )A(s ). (2.6) 
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The parametrization ( |2.5| ) allows one to write Eqs. fl2.1p - Q2.3| ) as a nonlinear, 
singular integral equation for the phase shift 5, 

——sin 25 s = a+ K - } - j --^L s £ 4, s ] , 

2a(s) n J 4 x—Ax—s 

u(x) = M*)-^*<») ( , 7) 



with boundary condition ) . We shall refer to Eqs. ( |2.1| )-( J2~1D or Eqs. 
and ( |2.7| ) as Roy equation with input (a, A). We assume in the following that 
this input is non-vanishing. 

Once a solution of the Roy equation is known, the real part of the amplitude 
above so is obtained from the dispersion relation (|2.1|) , and / is then defined on 
[4,oo). 

The above formulation of the Roy equation is used in the following section. 
There exists an equivalent approach, based on the just mentioned fact that the 
dispersion relation fl2.1|) can be extended to the half axis [4, oo). This implies 
that the amplitude / is the boundary value of an analytic function, holomorphic 
in C\[4, oo). The Roy equation then amounts to the construction of this analytic 
function. This method is described in detail in section 4 and used in sections 5 
and 6. 

The value of the phase shift at the matching point plays a crucial role in 
the following analysis || |K], 12|. However, the input absorptive part fixes S(so) 



through the boundary condition (|2.6|) only modulo it and up to its sign. In 
Fig. 0, we display sin 2 5(s ) as a function of the phase shift 5(s ), together 
with the quantity o~(sq)A(sq), shown with a horizontal line. The filled circles 
correspond to values of S(s ) that fulfill the condition ( |2.6| ). For reasons outlined 
PHI , we stick to phase shifts that are not too negative. Furthermore, the case 



in 



where the phase shift at the matching point is a multiple of tt/2 requires special 
considerations. This complication may be avoided with an appropriate choice 
of the matching point s in actual calculations. Unless stated otherwise, we 
therefore assume in the following that 

<K*o) > ~ , 

S(s ) ^ nvr/2 , n = 0,l,2,... . (2.8) 

The solution of the Roy equation is not unique in general ||, 0, |12| , and it is 
useful to divide the manifold of solutions into classes C„ , parametrized by the 
value of the phase shift at the matching point, 

7T 

: 5(s ) e {mT,mr + -) , 
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sin 2 5(s ) 




G(s )A(s ) 



271 



Figure 1: Solutions to Eq. Q2.6Q . The filled circles correspond to values of the 
phase shift S(s ) that fulfill the condition (|2.6| ). 



c; 



71 



5(s ) E (nn — — , nir) 



n 



0,1,2, 



(2.9) 



These classes are indicated in figure 0. To quote an example, a solution / with 



S(s ) 



5tt 



belongs to the class C± . 



Aside from the question of uniqueness, we will be concerned with the sin- 
gularity structure of the solutions: every solution of (p. 71 ) with arbitrary input 
(a, A) is regular on (4, sq), but singular and only Holder continuous at the end 
points of that interval. Whereas the singularity at s = 4 is due to the threshold 
behaviour, the one at the matching point is unphysical, because the position of 
So is arbitrary. It is only for the special class of analytic inputs - which will be 
specified in section 6 - that there exists a solution that is regular at Sq. 

Although the original system of coupled, nonlinear and singular integral 
equations has been reduced to the simplified equation Q2.7|), we are unable 
to solve it explicitly for an arbitrary input. For this reason, we have to rely 
on alternative tools to achieve our goal. In particular, we i) investigate the 
infinitesimal neighborhood of a given solution, ii) construct explicit solutions in 
the case where the matching point is moved to infinity, and iii) investigate the 
manifold of solutions for a conveniently chosen specific input. A clear picture 
of the multiplicity of the solutions will emerge in this way, and the role of the 
unphysical singularity at the matching point can be investigated in a satisfactory 
manner. 
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C + 



Class 



m 



5(s ) 







71 



2tc 



Figure 2: The classes as defined in ( |2.9| ). The number m in the lower row 
equals the dimension of the infinitesimal neighborhood of a given solution, see 
section 3. 



3 Linearization and existence of ambiguities 

We assume in this section that the Roy equation with input (a, A) does have a 
solution 5, and seek for solutions 8' that are nearby. These can be determined 
by linearizing and explicitly solving the integral equation for the difference of 
the phase shifts. As already mentioned in the introduction, this method has 
two apparent drawbacks: first, one does not prove the very existence of the 
solution 5. Second, it is not shown that the m-dimensional neighborhood of a 
solution is embedded in an m-dimensional manifold. On the other hand, as we 
will show in sections 4 and 5, we can construct explicit pairs (a, A) for which 
a solution is known and for which the ambiguities found below are present. In 
any case we believe that, despite its shortcomings, the linearization method is 
very useful and enlightening, in particular so in view of its simplicity. 

The solutions of the linearized integral equation show that, if a class with 
n > is non-empty, it contains a continuous family of solutions @. This 
multiplicity structure is identical to the one indicated in the work of Atkinson 
and Warnock fllP) , and the investigations of the nonlinear case carried out in 



later sections support this picture. 

We now show how this result is obtained when using the method of || [12] . 
By assumption, both phase shifts 5 and 5' satisfy the integral equation ( |2.71 ). 
We wish to determine the difference 

A(s) = 8'(s)-5(s). (3.1) 

Equation ( j2.7| ) for S and 6' results in a nonlinear singular integral equation for 
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A. Assuming A to be small, linearization of this integral equation gives 

,. c/ , u/ , , ,J p, 1 sin(25(x))h(x) [A . .„ . 

cos( 25 (s))h(s) = (s - 4)- Jjr dx - — v " K ' , s e 4, s , (3.2) 

71 J 4 X — 4 X — S 

where 

h( 8 ) = ^4- (3.3) 



To be consistent with the condition (|2.6|) , we require furthermore that 5'(sq) = 
6(s ), or 

h{s ) = . (3.4) 

The original Roy equation is replaced by the singular linear integral equa- 
tion (|3.2|), to be solved with the boundary condition ( |3.4| ). The latter shows 
that we can determine in this manner only those /' which belong to the same 
class as /. 



Constructing the solution of (|3.2| ) is equivalent to solving a boundary value 
problem for analytic functions - a so called Hilbert problem |[27| . To make 
this article self-contained, we discuss the procedure in appendix A, where it is 
shown that the general solution of (|3.2|) is given by M 



h(s) = (s - 4)G(s)P(s) 



with 



G(s) 



1 



exp 



s E [4, s ], 
5{x) ~ 



— £ dx 

IT J4 X — S 



(3.5) 
(3.6) 



and where P(s) is an arbitrary real polynomial of degree m — 1, with 

if - f < 8(s ) < |. 
The symbol [x] in ( p.7|) denotes the greatest integer not exceeding x. Polyno- 



m 



(3.7) 



mials of negative degree are considered to be identically zero here and in the 
following. For 5(sq) > , the number m coincides with the so-called index 
of the Hilbert problem. For a monotonically increasing phase, it counts the 
number of times 8{s) goes through multiples of tt/2 as s varies from threshold 
to the matching point sq- We indicate in figure [2] some of its values. 

The following proposition summarizes the results obtained in this section 
and in appendix A. 

Proposition 1 Let 5 be a solution of Eq. r \ ). It is an isolated solution of 
that equation if —\ < S(s ) < |. If 5(s ) > | ; the infinitesimal neighborhood 
of 6 is an m-parameter family of solutions 5' with <5'(s ) = S(s ), where m is 
given in (37i). 
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Comments 



An obvious question concerns the interpretation of the m parameters in 
the polynomial P in ( |3 . 5| ) . In the case where the phase shift 5 is monotoni- 
cally increasing and where m is an even integer, / exhibits m/2 resonances 
on (4, So). The nearby /' has also m/2 resonances, with slightly modified 
positions and widths. These changes are fixed by the m coefficients in 
P. This is verified in an example that we study in subsection 5.1. The 
situation is not so simple if m is odd. The case m = 1 is illustrated in 
subsection 5.3. 



Atkinson and Warnock [1(J have studied these parameters in a single 



channel Roy equation that has a structure which is similar to the one 
studied here. Using an iV/D-method, they find that, when m is an even 
integer, the parameters are related to the position and residues of the 
CDD poles between threshold and the matching point. In the case where 
m is an odd integer, one of these parameters is connected with the singular 
nature of the N equation. 

3. The difference A is singular at the matching point. This is a signal of 
the singularity of the general solution of the Roy equation mentioned in 
section 2. 

4. Our discussion deals exclusively with restrictions to the low-energy inter- 
val [4, So] • Once A is determined on that interval, the value of /' for all 
real s is obtained from the dispersion relation ( J2.1| ). In the linear regime, 

f W=n . ) + <iz£to a r ib *EVMm 1 . (3.8) 
7r e\o J4 x — 4 x — s — ie 

Consistency with unitarity is then by no means ensured above So: if erf 
stays on or inside the Argand circle, af may well be outside. This shows 
that physical requirements which are not encoded in the Roy equation can 
reduce the ambiguities. 

5. Our method is also applicable if the first inelastic threshold is below 
So, provided that the absorption parameter r\ is known. The represen- 
tation fl2.5Q is replaced by 

f(s) = r^pr [r](s)e^ - l] . (3.9) 

The form of the linearized equation ( |3.2| ) is unchanged, whereas rj enters 
into the definition of the unknown, h(s) = rj(s)a(s)~ l A(s). With this 
modification, (|3.5|) remains valid. 
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4 The full amplitudes 



4.1 Matching at infinity 



It is instructive to consider the case where the matching point Sq is moved to 
infinity, because the Roy equation can then be solved explicitly. The amplitude 
vanishes at infinity |L7| in this result of which the input reduces to 

the scattering length a. Equation (|2.7|) becomes^ 



2<t(s 



sin(2<5(s)) 



a + 



(s — 4) r°° dx sin 2 5(x) 



7T J4 x — 4 a{x){x — s) 



s e [4, oo] 



(4.1) 



to be solved for given scattering length a. For a > 0, a solution is provided by 
the phase shift of the amplitude 



fi(s) 



--Pis) 



-1 



where p is the Chew-Mandelstam function 



1 



p(s) = -a(s) Un 



a(s) 



71 



1 + <t(s) 



171 



s> 4 . 



(4.2) 



(4.3) 



In figure ^ we display the phase shift Si of fi with a solid line for a = 0.5. This 
is not the only solution - further examples are e.g. the phase shifts of 



1 
1 



s — 4 r 2 



fi 
1 

T + ~ 

J2 S 3 



S2 
S - 



-As- s 2 
4 r 3 



4s - s 3 



Si > 4 , n > . 



(4.4) 



The corresponding phase shifts are again displayed in figure |3], for (s 2 ,r 2 ) = 
(30, 40); (s 3 , r 3 ) = (45, 40). We note that the phase shifts tend to multiples of tt 
at infinity, 5j(oo) = (z — l)vr,i = 1,2,3. [The complications with S(s ) = rnr/2, 
that we mentioned in section 2, disappear when s = oo.] 

These examples show that the Roy equation with matching point at infinity 
allows for many solutions. The poles at Sj are CDD poles [p8fl , and one has to 
specify the phase shift at infinity as well as the CDD parameters T{, Sj in order 
to pin down a solution completely. The unsubtracted version of equation ( (HQ 
is discussed in detail in UTi 



2 A subtraction is in fact not needed here JlTj. We stick to the present formulation for an 
easier comparison with the Roy equation at a finite matching point. 
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Figure 3: The phase shifts #1,2,3 that correspond to the solutions /i,2,3 of 
equation ( [4.1| ), with parameters a = 0.5, (s2,r 2 ) = (30,40), (s3,r 3 ) = (45,40). 



4.2 Matching at finite energy 

We now present an approach that will allow us to construct exact solutions of 
the Roy equation with a finite matching point Sq. Although we cannot solve 
the Roy equation with an arbitrary input, we can construct explicit amplitudes 
/ that satisfy the dispersion relation ( |2.1|) and verify elastic unitarity below 
some Sq. An amplitude with this property defines an input cif = /(4), Af(s) = 
lmf(s),s > so, and is itself a solution of the Roy equation with this input. 
This holds true e.g. for the amplitudes fa in Eqs. (|4.2|) and ( fill) . Our goal is to 
show that - given / - one can construct other solutions /' ^ f, with the same 
input. In this manner, we find that the Roy equation has in general solutions 
in different classes C%. This cannot be seen in the linearization framework 
discussed in the last section, since there, one has to assume that the phase 
shifts of the original and of the new solution coincide at the matching point, as 
a result of which the old and new solutions stay in the same class. 

Let us describe the procedure in detail. Until now we worked with ampli- 
tudes that are complex functions of the real variable s > 4. As announced in 
section 2, we now define our amplitudes as analytic functions in the complex s- 
plane, cut along the real axis for s > 4. The corresponding amplitudes verifying 
Eqs. ( |2.1|) - (|2.4| ) are the boundary values / + of /, defined as 

/+(*) = Hm/(s + ie) , se[4,oo). (4.5) 
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In particular, we consider the set of functions / with the following properties: 



i) / is holomorphic in C\[4, oo) and verifies the dispersion relation (|2 
written for / + . 

ii) /+ is elastic below the matching point, 

f + (s) = -±-e^ sin5( S ) , s £ [4, s ] . (4.6) 

iii) /+ satisfies the regularity requirements listed in subsection A. 2. 

Let / satisfy i) — iii) - Its boundary value /+ is a solution of the Roy equation 
with input 

a f = 7(4) , 

A f (s) = Im/ + (s) , s>s . (4.7) 

In this and the following section, we show how to construct functions f ^ f 
that satisfy i) — iii) with 

/'(4) = a f , (4.8a) 

lmf' + (s)=A f (s) , s>s . (4.8b) 

It is clear that f' + is then also a solution of the Roy equation with input 
( |4.7|) , and the existence of an /' with the above mentioned properties therefore 
establishes the non-uniqueness of the solution of the Roy equation. 

In order to construct an /', we use the following Ansatz, 

1 1 +(.-4)S!> (4.9) 



/'(») f(s) v 'D(s 
Here, H is an Omnes-type function Q], an analogue of G defined in ( |A.9| ) 



H(s) = ( — - — ) exp 
'-so 



2 r°° dx 9(x) 



7T Jsq X X — S 



(4.10) 



with m given in (|3.7p . The function 9 is Holder continuous and equal to arg/ + 
modulo 7T, with 9(s ) = 5(s ). (We assume here in addition that Im / + > 0. As 
a result of this, we can define arg/ + such that < arg/ + < ir.) The function 
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D is meromorphicQ in C\[so, oo) - which ensures that ii) is fulfilled - and has to 
be constructed such that also the remaining conditions are satisfied. We write 

D = D 1 + D 2 , (4.11) 

where D\ is a meromorphic component of D, and where D 2 is regular in 
C\[sq, oo). Condition ( |4.8a| ) requires D(4) ^ 0. The definition of H has been 
chosen such that the condition ( |4.8b|) amounts to a simple linear constraint on 
D 



2; 



Im D 2+ (s) = fi{s) , s > s , 

K s ) = (s-A)\H(s)\A f (s). (4.12) 

Once this condition is fulfilled and D(sq) is not zero, Im /' is continuous at 
s if H(sq) = 0. This holds true if S(s ) > 0. Choosing a particular function 
D 2 verifying the condition ( 4.12|) , the arbitrariness in /' is entirely contained 



in D\. This function has to be such that /' satisfies i) and iii). Therefore, 
1/ /' and have to be nonzero on C\[4, oo) and on [4, oo), respectively. The 
non-uniqueness of the solution of the Roy equation is due to the very existence 
of such functions D\ - we provide explicit examples in the following section. 
Our discussion leads to 

Proposition 2 Let f be an amplitude verifying conditions i)-iii), with Im/ + > 
0, and 5{sq) > 0. Let f ^ f also verify i)-iii), together with the conditions 
(\4-8j) . Then f can be written in the form - ( \4-H\ )> where D\ is a mero- 



morphic component of D, and D 2 is regular in C\[s , oo), with spectral function 

In the N/ D approach , an important part of the non-uniqueness is due to 
arbitrariness in the CDD poles. Their relation with the function D in Eq. ( [4.11 
is explained in appendix B. 



5 Explicit exact solutions 



Here, we illustrate the procedure described in subsection 4.2 with specific ex- 
amples. An amplitude fulfilling the conditions i)-iii) of that subsection can be 
written as follows, 



i r 1 

- + (s - 4 0s - p(s) 
a 



(5.1) 



3 In order to simplify the presentation, we shall - without any further mention - use the 
fact that all our analytic functions satisfy F(s) = F(s). 
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where is a suitable function that is meromorphic in the complex s-plane, cut 
along the real axis for s > sq. The amplitudes fi in Eqs. ( |4.2|) and (|4.4j) have 
this form with rational 0. To keep our calculations simple, we work with the 
resonant amplitude f 2 in (|4.4|). We drop indices and write 

-l 



1 s — 4 r . . 

- + — ^_ p(s) 



(5.2) 



where s p > 4, r > 0, a > 0. As / is elastic on [4, oo), the argument 6 in ( |4.10|) 
coincides with the phase shift 5. In figures ^ and [5], we display the behaviour of 
/ in the case where 

a = 0.5 , s p = 40, r = 25 . (5.3) 

The imaginary part of af is shown in Fig. |], whereas its phase shift 5 is indicated 
with a solid line in Fig. Note that the matching point sq does not occur in / 
- it may be chosen at one's convenience. We observe that s p becomes a CDD 
pole - as defined in |I| - if s p < sq. 



In the following, we keep the amplitude / fixed and bring it into class 
or Cf by an appropriate choice of sq. To illustrate, for so = 50, / is in class 
Ci , and s p is a CDD pole. 

We notice that / is a solution of the Roy equation with input (a/, Af) which 
is regular at sq. We are thus dealing with one of the special inputs mentioned 
in section 2. 



5.1 Shift and suppression of a resonance: / G Ci, f G 

°i ? °o 

We place s above s p , as a result of which n < 8(s ) < 3n/2, such that / 
belongs to class Ci . We first construct amplitudes /' in the same class Cf. 
This implies that S'(s ) = S(s ). We shall find that, qualitatively, /' is obtained 
from / by a shift of the position and a change of the width of its resonance. As 
S'(sq) G (it, 37r/2), 1/f has a single pole on (4, s ). The second term in the 
Ansatz ( J4.9| ) has to cancel the pole of 1// at s p and replace it by a new pole at 
ftp. 

It is convenient to redefine D\ and D 2 in ( [4.1 1|) by writing 

D(s) = (s- s p )(s - s'^D.is) + D 2 (s)} , (5.4) 
where D\ is meromorphic, and 

D 2 (s) = - r dx- ^ ■ (5-5) 

7T J so [X — S p ) [X — S' J X — S 
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Figure 4: The imaginary part of the function af in (|5.2| ). The parameters 
used are the ones in Q5.3|). The phase shift 5 of / is displayed with a solid line 
in Fig. m 



The integral converges because fj,(x) = O(x) at infinity. We require regularity 
of 1/ /' at s p and fix the residue r' of its pole at s' p , with 

Sp ^ s p , 4 < s p , s' p < s , r' > . (5.6) 

This gives two conditions which completely determine a two-parameter Ansatz 
for D\. One finds that the adequate Ansatz is 

D x {s) = . (5.7) 

as + p 

The two constraints on D\ give 

<*=---, P= f -i-^, (5.8) 
R' R R R' y 1 



with 



_ 
r 



(sp-4)H(s P Y (s' p -4)H(s' p y (59) 

R=l + f(s p - s' p )D 2 (s p ), R' = l + f'(s p - s' p )D 2 {s'p) . 
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The function H is obtained from (|4.10 ) with m = 2. Notice that a and (5 are 



small if the pair (s' p ,r') is close to (s p ,r). In particular, let 

e = max (\izM } ij^zil) f (5.10) 



with e small. We then have 

a = 0(e), /3 = 0(e). (5.11) 
Inserting the expressions (|5.4j) and ( |5.7|) into equation (|4.9| ) we get 



1 1 | (s-4)ff(s) aa + g ^ 



/'(s) /0) (s - s p )(s - 1 + (as + P)D 2 (s 

We see that /' is regular in C\[4, oo) and is bounded on its cut if 1/ /' has no 
zero. Eqs. (|5.12 ) and ( 5.11j) show that 1/f is close to 1// if e is small and s is 



outside the vicinity of s p and s' and not too large. The inverse amplitude 1/f 
is nonzero for such values of s because 1/f is nonzero. One finds that 1/f is 
also nonzero near s p and s' p . Since D 2 and H behave at infinity as s _1 In |s| and 
ln 2 |s|, respectively, the second term in (|5.12|) is no longer a small correction 
when \s\ becomes large, and a detailed analysis is needed. It is relatively easy 
to see that no unwanted zeros show up if e is small and (3 /a is not too large. 
Consequently we are sure that under these conditions ( |5.12| ) provides a two- 
parameter family of solutions of the Roy equation. This example provides an 
illustration of the N/D framework outlined in appendix B. In that language, s p 
and s' p correspond to zeros zj. The CDD pole of / at s p is removed and replaced 
by a new one at s' p . 

We have numerically verified that e need not be small for the above conclu- 
sions to hold. In particular, there is a finite interval for e such that i) there are 
no additional singularities in /', and ii) the dispersion relation (|2.1| ) is fulfilled. 
To illustrate, we choose 

s' p = 35 , r' = 20 , s = 50 . (5.13) 

In this case one has, using for / the parameters displayed in ( |5.3|) . 

\s'„ — s v \ \r' — r\ 

0.125 , l - = 0.2 . (5.14) 



s p r 



These quantities are not small. Nevertheless, /' is a solution of the Roy equation 
in class C^. We display the phase shifts of / and /' in Fig. || with a solid and 
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Figure 5: A shift in the class C± . The solid line corresponds to the phase shift 
5 of the function / shown in Fig. f|, whereas the dot-dashed line displays the 
phase 5' of /', evaluated from ( |5.12| ), with a, (3 calculated from ( |5.8| , |5.9| ), using 
the parameters Q5.13 ). 



dot-dashed line, respectively. As required, the two phase shifts agree at the 
matching point s - the parameters ( |5.13| ) really correspond to a shift Ci — > C^. 

An expression for /' can be obtained in the linear regime from the solution 
fl3.5p of the problem on [4, s ]. It coincides with the result ( |5.12j ) to first order 
in e, if the polynomial in (\3.l 



is 



P(s) = C[(f - f)s + fs' p - f's p ] 



(5.15) 



where C is a constant determined by the phase shift 5. This establishes the 
relation with the coefficients in the arbitrary polynomial P in our example and 
provides an interpretation of these coefficients. 

The method allows the construction of solutions /' which are no longer res- 
onant. This is simply achieved by setting the residue r' to zero. Equation 
then gives as + (3 = —f(s — s' p )/R, and 1/ /' becomes 



1 1 _s-A 

— -, — = h T 

f'(s) a s - s p 



His, 



H{s} 
Lis) 



(5.16) 
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where 



> - s p )D 2 (s) + (s p - s' p ) (D 2 (s) - D 2 (s p ))\ . (5.17) 

As 1/f has no pole on [4, sq], S'(so) = <5(so) — ir, and /' G Cq. We have to 
make sure that 1/f has no unwanted zeros. The conclusions reached in the 
case r' 7^ apply here and 1/f is non-vanishing if the residue r is small, i.e. if 
the resonance in / is narrow. On also finds that Re(l//') does not have a 
zero on [4, s ] either. Therefore the phase shift 5' stays below ir/2, and /' is 
non-resonant. It is easy to see that L does not depend on the choice of s' p . 
This means that the amplitude /' defined in ( |5.16|) effectively contains no free 
parameter. This amplitude is the unique solution of the Roy equation with input 
( f4.8|) belonging to class Cq . This uniqueness is in accordance with the results 
of section 3. In the N/D language, the uniqueness stems from the fact that /' 
has no CDD pole, and 5'(sq) < ir/2. What is new with respect to section 3 is 
that our example shows explicitly that the same Roy equation has solutions in 
different classes C^. At given input (a, A), the number of resonances below s 
is not fixed unless one imposes the precise value of S'(s ). 



5.2 Implantation of a resonance: / G C\ , f G C2 

We corroborate our last statement by examining another Ansatz for Di in the 
context of the last subsection. We replace (|5.7|) by 



D x {s) = -(s - si) , (5.18) 
a 

and again require that 1/f obtained from Eqs. Q4.9D and (|5.4j) be regular at s p 
and have a pole at s' p with residue r', with conditions ( p.6|) . One finds 

1 — If I \2 

si = ^(s' p + s p )-^(s' p -s p ){f' + f-f'f(D2(s' p ) + D 2 (s p ))(s p -s p )} , 
AT = r -f' + ff'{s p -s' p ){D 2 {s' p )-D 2 {s p )} . 

(5.19) 

We see that in general, a = 0(e), whereas si is anywhere on the real axis. 
Equation Q5.12|) is replaced by 

1 _ 1 , (s-A)H(s) a 



f(s) f(s) (s - s p ) (s - s' p ) s - si + aD 2 (s) ' 
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If e is small, 1/ /' has a pole s"i near s\ with small residue, 0(e). This implies 
that 1/f exhibits a zero near s\ which can preclude the Ansatz ( |5.18| ). The 
discussion becomes delicate if s± is close to 4 or so and the next statements are 
valid if Si is outside 0(e) neighborhoods of 4 and s : 

i) if s± < 4 or s\ > sq, f has a pole on the real axis near si; 

ii) if 4 < si < so, /' has a pair of complex conjugate poles close to si. These 
poles are in the first sheet of the branch point s = 4 if a < 0; they are in 
the second sheet if a > 0. 

We see that the Ansatz ( |5.18| ) has to be rejected for pairs (s' p , r') such that s"i < 4 
or &i > so. It is also inadequate if 4 < si < s and a < 0. These requirements 
on a and s± are fulfilled in a wedge-shaped domain of the (s' p , r')-plane, with 
apex at (s p , r). 

If a > 0, /' generates a solution of the Roy equation in class C£ (5'(sq) = 
S(sq) + 7r). This solution displays two resonances on [4, s ] : a shifted resonance, 
(r, s p ) — > (r', s' p ), and an implanted narrow resonance near s±. We illustrate this 
transition in figure |6]. There, we display the phase shifts of / and /' with a solid 
and a dot-dashed line, respectively, for 

s ' p = 47 , r' = 20 , s = 50 . (5.21) 

These parameters are again outside the linear regime. The difference of the 
phase shifts is exactly it at the matching point So = 50. This amounts to a 
special case where the CDD pole at s p is replaced by two poles at s' p and s"x. 

We have pointed out in section 3 that the new solutions exhibit a singular 
behaviour at the matching point sq. We illustrate this feature in figure ^|, where 
we display the real and the imaginary part of / and of /' with a solid and a dot- 
dashed line, respectively, for the situation displayed in figure |[ The singular 
behaviour of /' is manifest, both in the real and in the imaginary part. Note 
that, since the imaginary parts of / and /' agree for s > so, the real parts also 
agree at s = so due to unitarity. Above the matching point, they are, however, 
in general different. Because the real part of / is positive in the present case, 
the figure shows that /' is inside the Argand diagram, at least in an interval 
above sq. 



5.3 Shift of a resonance: /, /' G C 1 



In this example we bring / in ( |5.2|) into C t by pushing sq below s p in such a 



manner that the resonance position is still below sq, i.e. ir/2 < 8{sq) < n. [The 
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Figure 6: The phase shifts corresponding to a shift — > C^. The solid line 
is the phase shift 5 of the function / shown in Fig. ^, whereas the dot-dashed 
line displays the phase shift 5' of /', evaluated from (|5.20 ). The parameters a 
and Si are evaluated from ( |5.19| ), with (|5.21 ). 
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Figure 7: The cusps in /'. We display the real and imaginary parts of / and 
/', for the situation displayed in Fig. ||. 
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pole at s p is no longer a CDD pole.] The function H used in the previous two 
subsections has now to be replaced by a new one that we call H. It is given by 
( [HOp , with m — 1. Correspondingly, ii defined in ( |4.12| ) becomes fx. It behaves 



as s at infinity and 3 subtractions are needed in the construction of D 2 . A 



convenient redefinition of D\ 2 in ( 4.11j) leads to 



where 



D{s) = {s-s p )[D 1 (s) + {s-s p )D 2 {s)] , 
dx p,(x) 1 



Do(s) 



s— 

7T j s 



X [X — S p )* X — S 



(5.22) 



(5.23) 



The integral converges because fx has a second-order zero at s p . As before, we 
require 1/ /' to be regular at s p . This implies 



f + 



H(s p ) 
D 1 (s p ) 







(5.24) 



with f = r/(s p — 4). This equation is consistent with real values of f and D\(s 
because H(s p ) is real. A suitable Ansatz for D± is 



DAs) 



H(s r 



r(s„ — s' 



(5.25) 



Condition (|5.24 ) is fulfilled with s' p a free parameter, s' p > so- Insertion of 
Eqs. flOp and (|5T25|) into Eq. gj) gives 



1 

a 



4)r 



1 - 



F(s) 



(5.26) 



where 



L(s) 



H(s p ) s-s> p + L(s) 



l+r(s' p - s p )D 2 ( 



(5.27) 



Using e defined in ( |5.10| ), we have L = 0(e 2 ) as long as s is not too large. For 
such values of s, 1/ /' is close to 1/ f outside a neighborhood of s p and s' . We 
notice that Im L + ^ for s > sq and 1/ /' has no real pole whose real residue r' 
could be imposed (there is a pair of complex conjugate poles near s' p located in 
higher sheets of the logarithmic branch point at sq). One finds that 1/ /' has no 
unwanted zeros for any s if e is small enough, and we end up with a family of 
solutions of the Roy equation in Cf indexed by the single parameter s' p . This 
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is in accordance with section 3. The shift of the resonance position and the 
change of its width are now correlated and fixed by the value of (s' p — s p ) . Both 
/ and /' are without CDD poles, and the non-uniqueness is due to the fact 
that 7r/2 < 8(s ) < 7T. 



6 How to ensure uniqueness 

The solution of the Roy equation is not unique - the behavior of a partial wave 
below So cannot be predicted in a unique way from an input (a, A). Information 
on the value of the phase shift at Sq is needed and free parameters have to be 
fixed if S(sq) > 7t/2. This non-uniqueness restricts severely the efficiency of a 
Roy equation as a tool for the construction of a low-energy extrapolation. We 
show in the following, how uniqueness can be restored in principle by imposing 
additional physical requirements. The unphysical singularity at the matching 
point is removed at the same time, as it has to be. 



6.1 Examples of unique solutions 

We start with three illustrative examples. 

Example 1 In section 5 we have constructed, starting from the amplitude 



(|5.2|) , several functions /' that satisfy the Roy equation with matching point 
so = 50 and input generated by /, 

a = /(4) ; A(s) = lm f(s) , s>s . (6.1) 

Two cases are displayed - in terms of their phase shifts 5' - in Figs. |5] and || As 
we already mentioned, these phase shifts develop singularities at the matching 
point, see Fig. [?|. Below we will show that, had we required the new solution 
/' to be regular at s , there would be exactly one solution of the Roy equation 



with input (dTT), namely / itself. In other words, all solutions f ^ f develop 
singularities at s . It turns out that this property of the input (|6.1|) is due to 
the fact that the amplitude is elastic above the matching point. Whereas / is 
elastic on the whole interval [4, oo), the following example shows that a finite 
interval containing [4, s ] suffices to render the solution unique. 

Example 2 We consider the function /' constructed in subsection 5.2, see 



Eqs. ( 5.20|) , (|5.21 ) and Fig. [| dot-dashed line. Suppose we wish to construct 
solutions of the Roy equation with displaced matching point s' = 45 and input 
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denned by the scattering length and the absorptive part of /', 

a = f'(4) ] A(s)=lmf\s),s>s' . (6.2) 

In the language of section 2, this problem belongs to class C^'- the phase shift 
at the matching point is 37r/2 < S'(s' ) < 2n, see Fig. ||, and according to 
proposition 1 in section 3, the Roy equation with input (|6.2|) has therefore a 
3-dimensional manifold of solutions. The uniqueness statement in example 1 
is also true here: Suppose we seek for solutions of the Roy equation with input 
(|6.2|) and require that the solution is regular at s = s' . There is again exactly 
one solution, namely/' itself. 

It is obvious that we are dealing with special inputs - we call them analytic 
inputs below. The following example displays an input that is not analytic. 

Example 3 Consider again the amplitude (|5.2|). We set the residue r to zero, 

a = /(4) + e ; A(s) = f(s) ,s>s ;r = 0, (6.3) 

with So = 50. For sufficiently small e ^ 0, one can show that there is no solution 
that is analytic at s , see subsection 6.3. 



6.2 Analytic input and uniqueness 

To generalize our findings, we exploit smoothness properties of the amplitudes 
which are not explicit in the Roy equations. The partial wave amplitudes en- 
ter this equation as boundary values of analytic functions. Boundary values 
need not be smooth, and this is compatible with the dispersion relation (|2.1| ). 
However, it is quite remarkable that smoothness is imposed by elastic unitarity: 



Proposition 3 Let f be regular in the complex s-plane, cut along the real axis 
for s > 4, and let its boundary value f + verify elastic unitarity on [4, s\] . Then 
the real and imaginary parts of f + are separately holomorphic in a complex 
neighborhood of (4, s\). 



A proof of the proposition is given in PBJ. Notice that the parameter s\ 



need not coincide with the first inelastic threshold Si ne i - the proposition is true 
for any si > 4. Taking si = s , one concludes that all solutions 5 of the Roy 
equation ( p.7| ) are regular in the interval (4, s ). 

Consider amplitudes that satisfy the conditions of proposition 3 and in addi- 
tion verify the Roy equation with input (af,Af(s)) = (/(4), Im/ + (s)) for some 
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so < si- The proposition tells us that / is regular at sq. A second solution of the 
same equation may be singular at sq, because it need not be elastic on [sq, si]. 
The following proposition shows that / is in fact the only solution which is 
regular at Sq. This is an important feature, allowing one in principle to identify 
the physical extrapolation within the manifold of solutions. 



Proposition 4 Let f be an amplitude that satisfies the conditions of proposition 
3 and that furthermore verifies the Roy equation with input (af,Af) for some 
So < S\. Let f'^fbea second solution with the same input. Then f is regular 
at sq, whereas f is singular. 



We relegate a proof of the proposition to appendix C. It elaborates obser- 



vations made in || and JTJj. The examples 1 and 2 given above fulfill the 
conditions of the proposition with si = oo and si = 50, respectively. Therefore, 
according to proposition 4, there is exactly one solution of the Roy equation 
that is regular at the matching point. 

Proposition 4 tells us that there is a special class of inputs that allow a 
unique solution that is regular at the matching point Sq, as announced in section 
2. According to proposition 3, the high-energy absorptive part of an input 
belonging to that class has an analytic continuation from [s , «i) into a complex 
neighborhood of (4, si). For this reason, we say that the members of our special 
class are analytic inputs. A physical amplitude / defines an analytic input 
(a/, Af) if so < Si ne i. Uniqueness is achieved in the sense that the corresponding 
Roy equation has exactly one solution, coinciding with /, which is regular at 

Sq. 

The existence of a class of inputs ensuring uniqueness and regularity at the 
matching point has been established in an indirect way. We have no direct and 
complete characterization of an analytic input. There are involved constraints 
apart from analyticity of the high-energy absorptive part A, and one cannot 
decide directly if a given input is an analytic one. In particular, the scattering 
length a is fixed by A and is not an independent parameter, see example 3 
above. We arrived at a unique solution of the Roy equation by choosing inputs 
which are compatible with elastic unitarity above the matching point. In the 
physical context, this means Sj ne ; > s . One can prove that uniqueness is also 
obtained if sq > Si ne i, provided that the inelasticity is sufficiently smooth. 
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6.3 Approximate input and cusps 



Although we are working here with model amplitudes without left-hand cut, 
the results of subsection 6.2 hold true in the physically realistic situation with 
left-hand cut. As a physically relevant input is analytic, with s\ = Sj ne ;, we 
conclude that it allows a unique solution regular at s if s < s ine i. However, a 
physical input is only approximately known, and one is faced in practice with an 
arbitrary input, as a result of which non-uniqueness and unphysical singularities 
do occur. This is the reason why we have analyzed in detail the Roy equation 
in its general setting. 

If we know that a high-energy absorptive part A belongs to an analytic 
input, but the corresponding scattering length is not precisely known, one has 
to work with a trial input (a', A), with a' 7^ a. We expect that all the resulting 
solutions will have cusps at sq. This can be established if a' is infinitesimally 
close to a by using the techniques developed in j9|, |I2]|. For instance, if S(so) < 
n/2 and S'(sq) = S(sq), the neighboring solution /' is unique and its phase shift 
is given by 

6 , (s) = 6(s) + a(s)^(a'-a), (6.4) 

with G defined in ( |3.6|) . The ratio (5' — S)/(a f — a) is shown in Fig. § for a 
neighboring solution of the amplitude / in (|5.2| ), using 



Sp = 40 , r = 100 , s = 23 . (6.5) 

The cusp at So is clearly visible and the effect of a modified scattering length 
extends over the whole interval (4, Sq). If an incorrect value a' of the scattering 
length also entails a cusp when it is not close to a, the correct scattering length is 
specified by the fact that it defines an input allowing a unique solution without 
a cusp: one solves the Roy equation with several trial a', and uses 

a = a 'nocusp ■ (6-6) 

In this sense the scattering length of a physically realistic amplitude can be 
predicted in principle. However, even if we are sure that the high-energy ab- 
sorptive part is analytic, we do not know it exactly in practice, and we have to 
use an approximate form A'. We have to work with an input (a', A') which is 
meant to approximate the analytic (a, A). It is unlikely that there will still be 
a value of a' which removes the cusp in one of the solutions of the resulting Roy 
equation. What one can try in practice is to find the value of a' which minimizes 
the size of the cusp. In the coupled channel case, an analogous procedure may 
be used to avoid solutions that generate a cusp @. 
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Figure 8: The quantity (6' — 8) j{p! — a) according to Eq. (|6.4j) . The reference 
amplitude / is the one in ( |5.2|) , and the parameters used are given in (|6.5). The 
cusp generated at s is very clearly seen. 



7 Summary and conclusions 

The following points summarize the content of this article. 

1. In view of forthcoming applications [|J of the Roy equations M in the 
analysis of decays, we have considered here the one-channel Roy equa- 
tion. We have analyzed the multiplicity and singularity structure of its 
solutions for a given input (a, A) of scattering length a and high-energy 
absorptive part A. 

2. First, we have investigated the infinitesimal neighborhood of a given so- 
lution 5 ||. According to proposition [I] in section 3, this neighborhood 
contains an m-parameter family of solutions, where 

*w > f, (7-1) 

- f < S(s„) < f . 

The symbol [x] denotes the greatest integer not exceeding x. For a mono- 
tonically increasing phase, m counts the number of times S(s) goes through 
multiples of it/ 2 as s varies from threshold to the matching point Sq. This 
result illustrates that a given input (a, A) does not, in general, fix the so- 
lution uniquely. One has in addition to fix the phase shift at the matching 
point, and to determine the corresponding m parameters by other means. 
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3. Using proposition 2 of section 4, we have constructed in section 5 - starting 
from a given solution / with input (a/, Af) - additional exact solutions 
f f with the same input. The function /' contains in general several 
arbitrary parameters that may be used to either change the position and 
residues of the poles present in the inverse amplitude 1//, or to remove 
(or implement new) poles. This illustrates that a given input allows for 
solutions with a different value of the phase shift at the matching point - 
these phase shifts only have to satisfy the boundary condition ( |2.6| ). 

4. The solutions f ^ / so constructed have the property that they contain a 
cusp in the real and imaginary part at the matching point. An example is 
displayed in figure |7|. The origin of these cusps is made clear in proposition 
| of section 6. 

5. In case we know that a given absorptive belongs to an analytic input (see 
section 6 for this notion), we expect on the basis of proposition [|that the 
corresponding scattering length can be determined as the one that results 
in a solution of the Roy equation without a cusp at sq. 

6. Propositions 1 and 3 in sections 3 and 6 were established long ago in 
Refs. and P6| , respectively. As far as we are aware, propositions 2 and 
4 in sections 4 and 6 are, on the other hand, new results. 

In conclusion, if we know that we are dealing with an analytic input, the 
corresponding amplitude is the unique regular solution of the Roy equation. 
Because we are in general forced to use an approximate input, non-uniqueness 
and unphysical singularities do show up. A visible cusp in a numerical solution 
of the Roy equation with < 5(sq) < vr/2 is a signal that the input is not 
physical or a poor approximation of a physical one. The deficiency may be 
hidden in the scattering length a, in the absorptive part A, or in both. 
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A Solving the integral equation ( |3T2l) 



We first specify the regularity properties of the partial waves and phase shifts 
considered in the text. Then we construct the general solution of the integral 
equation 



A.l Holder continuity 

The class of Holder-continuous functions is the appropriate space to consider, 



because Cauchy-integrals map this space essentially into itself |27| . We explain 
this notion. 

Consider a complex-valued function / of a real variable x. The function / is 
called Holder-continuous in the interval [a, b] with exponent //, where < \i < 1, 
if there is a constant C such that 

1/0*0 -/(j/)|<C|z-y|" ; x,ye[a,b]. (A.l) 

We call these functions if -continuous for short, and denote the corresponding 
space with H% b . 

Proposition A Let f E H% b , with /(a) = fib) = and /i < 1. Then the 
function 

g{y) = fdx^- (A.2) 
J a x — y 



is also an element of H^ b , with the same exponent //. The same is true for 

g±(y) = lim t dx f ^ . (A.3) 



The proof may be found in [27 . 



A.2 Regularity requirements 



The input absorptive part A in (|2.1|) is assumed to be bounded and Holder con- 



tinuous in any finite interval [so, s' ] above the matching point sq. Furthermore, 
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we seek solutions of the Roy equation that are if-continuous in [4, sq] and have 
normal threshold behaviour (which implies /x < 1/2): 

i) A G ., , s' n > s ; A is bounded 



ii) / G H? 



4,so 

2 



iii) / = a + ia 2 a(s) + 0(s -4) , s \ 4 . (A.4) 



The phase shifts can then be chosen if-continuous as well, 

iv) 6 G H^ SQ 

v) 6 = a(s)[a + 0(s - 4)] , s \ 4 . (A.5) 



A. 3 Solution of the integral equation (|5.2| ) 

We first show that solving the integral equation (|3.2|) with the boundary con- 
dition fl3.4| ) is equivalent [E7 to solving a boundary value problem known as 
Hilbert problem. We then construct the general solution of the latter. 

We start by introducing an auxiliary function $ || [T2"|] , 

1 /" s o sin(25(x))/i(x) 



,1 /" s o sin rzdix ))ti (x) 1 /* \ 

z-4)-/ dx — v v v ; _ Ag 

^ 71 J 4 X-4 £ - Z 



For convenience, we use in this subsection the variable z to indicate complex 
values of the variable s. The function $ has the following properties if h is a 



solution of Q and fl3j): 



i) $(2) = $(2) is regular in C\[4, sq] . 

ii) The boundary values $±(s) = lim$(s ± ie) are .ff-continuous in [4, s ] . 

iii) $+(s) = e 4i5 ^$_(s) ; s G [4, s ] . 

iv) $(4) = . 

v) $(s ) = . 

vi) $ is bounded at infinity . 

Out of these, we only prove property ii) - the remaining ones are easy to verify. 
The unknown h in (|A.6| ) is an element of H% so , with < \l < 1/2. This 



follows from Eq. ( |A.5| ) - the same equation shows that h = 0(s — 4) at s = 4. 
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Together with h(so) = and with proposition A in subsection A.l, it follows 
that $±(s) G H£ SQ , with the same exponent \x. 

The problem to determine functions &{z) with i)-vi) is called a (homoge- 
neous) Hilbert problem. We conclude that each solution of the original integral 
equation generates via (|A.6| ) a solution of i)-vi). Vice versa, each solution of 
the Hilbert problem has a representation of the form ( |A.6| ), where h is given by 



h(s) = e- 2i5(s) $+(s) . (A.7) 



Furthermore, the function h is real and solves the integral equation (|3.2| ). To 
prove this, we write a Cauchy representation for — $(4)]/(z — 4), with a 

path that wraps around the cut [4, s ] . We then deform the outer part of the 
path towards infinity, in such a manner that only the integral above and below 
the cut [4, So] survives. By use of condition iv), we find 

$ (z) = (i_z3 r° dx (A8) 

1 ; 2ni J a x-4 x-z ' y ' ' 

Using properties i) and iii), the claim is easily proven, and we conclude that 



solving (|3.2| ) with (|3.4| ) indeed is equivalent to solving the Hilbert problem i)-vi). 



It remains to construct the general solution of i)-vi). First, we observe that 
the Omnes-type function M 



G(z) = ^— exp 

(s - z) m 



s o , S(x) 
dx 



(A.9) 



7T J4 X — Z 

with m defined in Q3.7|) , satisfies property i). The behaviour at sq is given by 
G(z) ~ ( So - z)\ 7 = -8(s ) - m , (A.10) 

— 1 < 7 < 1. The if-continuity and threshold behavior of 5 imply that G 
satisfies ii) and iii) except, possibly, at so- Outside so, (so — z) m G and its 
boundary values are nonzero and the function 

f w=( s ° ( A - n ) 

is regular in C\{so}. For z ^ so it is given by its Laurent series. The principal 
part of this series is identically zero because condition v) and Eq. (|A.10|) imply 
that F{sq) = 0. F(z) is therefore an entire function. Condition vi) tells us that 
yTn+ at infinity, result of which 

F(z) = (s -z)Q(z), (A.12) 
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where Q(z) is a polynomial of degree m. Condition iv) imposes Q(4) = 0, and 
this gives Q = if m = 0: the conditions i)— vi) allow only for the trivial solution 
in this case. If m > 0, Q(z) = (z — 4)P(z), with P a polynomial of degree m — 1. 
We conclude that the general solution of i)-vi) is given by 

m— 1 

<S>(z) = {z- A)G(z) c n z n , c n e IR , (A.13) 

n=0 

with $ = for m — 0. Using ( |A.7| ), the result ( |3.5| ) follows. 



B Connection with the N/D approach 



The N/D method |H| transforms the nonlinear Roy equation into a linear N 



equation, whereas the method displayed in section 4 linearizes the construction 
of the solution /', once a first solution / is known. In this appendix, we establish 
the relationship between the two methods. 

We write the N/ D representation of / as follows, 

<■»> 

The A^-function n is holomorphic in C\[so, oo), and the D-f unction d is holo- 
morphic in C\[4, sq], with possible CDD poles on the cut. Similarly, /' = n'/d'. 
We now determine the relation between the pairs (n, d) and (n', d'). 



Equation (|4.9Q gives 

m = y^m , (b.2) 

where 

F(s) = D(s) + (s-4)H(s)f(s) . (B.3) 

Whereas D is meromorphic in <C\[so,oo), F is meromorphic in C\[4, Sq], This 
crucial point is a consequence of Eqs. ( |4.10|) -( [4.12j ), which imply that ImF + (s) 
= for s > Sq. In order to turn flB.2j ) into a N/D representation, we need the 
set {pi} of poles of D, as well as the set {zj} of its zeros on [4, sq], and define 

D{s) = Y[{8-p i )Y[ 7 ±—D(8). (B.4) 
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To simplify the argument, we assume in the following that the iV-functions n 
and n' have no zeros on [4, so] - the discussion is more involved and the defini- 
tions (|B.5|) have to be modified if this is not the case. The N/D representation 
of /' is then obtained by eliminating D in favour of D in the expression 
and writing 



n 



(s) = Cn{s)D(s) ; d\s) = Cd(s)F(s) 



where C is a constant and where 



F(s) = D(s) + (s- A)H(s)f(s) I] — - II(« - Pi) ■ 



j S ~ Z 3 i 



(B.5) 



(B.6) 



The functions n' and d! defined in (|B.5|) have the correct analyticity properties. 
Denominator functions are normalized to one at infinity jnj. We have checked 
in our examples that F has a finite, nonzero limit at infinity. Therefore, d' is 
properly normalized if C — l/F(oo). For S(sq) > 0, d can be written as 



d( S ) = n 



k=l 



S 



s - s k 



exp 



, 8(x) 
dx 

71 J i X — S 



(B.7) 



where the Sk are the r CDD poles of /, s k G (4, s ),r = [5(s )/vr]. We see that 
d' is finite at Sk if Sk coincides with a zero Zj of D, and if F(zj) = 0. Otherwise, 
every s k and Zj IS 3j CDD pole of /'. Within our assumptions on n and n', the 
fate of the CDD poles is dictated by the zeros of D: these produce new CDD 
poles in /' or remove CDD poles present in /. The poles of D are points where 
/' = /•' 



C Proof of proposition 

Using the framework of section 4, we show that /' coincides with /, if Re f' + 
is regular at so [in the sense that it has a holomorphic extension into a circle 
of radius e and center s ] . See figure |] for the analyticity domains used in the 
proof. 



i) The dispersion relation (2.1) written for /' determines Re f, on IR. In- 



version of this relation gives 

Im f'Js) = -(s-A)-idx^- , (C.l) 

7T JIR X — S 

where <j>(s) = (Re f' + (s) - f' + (4))/(s - 4). The holomorphy of Re f' + in C 
implies the holomorphy of Im f' + in the same circle. 
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Figure 9: Analyticity domains referred to in the proof of proposition 4. 



ii) According to proposition ^, Im f' + has a holomorphic extension into a 
neighborhood N' of (4, sq), with so on its boundary. Combining this with 
the previous result we see that Im f' + is holomorphic in N = N' U C, a 
domain extending up to so + e. 

iii) Proposition 3 tells us that the high-energy absorptive part A, originally 
defined on [so, oo), has an analytic continuation into a neighborhood N 
of (4, si), which coincides with Im/ + on [4, s ]. As /' is a solution of the 
Roy equation with input (a, A), we have 

Im f' + (s) = A(s) (C.2) 

for s < s < s + e. In view of this equality and of the regularity of Im f' + 
in N, Im f' + has to be equal to A on [4, s ]. As A is equal to Im / + on 
that interval, we conclude that /' = /. 

iv) Similarly, if Im f, is assumed to be regular at Sq, one concludes that Re / + 
is regular at So and /' = /. 

Therefore, /' has to be singular at s if /' ^ f. This is the content of 
proposition 
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